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Abstract 

The propagation of intense acoustic waves in a one-dimensional phononic crystal is studied. The medium consists in 
a structured fluid, formed by a periodic array of fluid layers with alternating linear acoustic properties and quadratic 
nonlinearity coefficient. The spacing between layers is of the order of the wavelength, therefore Bragg effects such 
as band-gaps appear. We show that the interplay between strong dispersion and nonlinearity leads to new scenarios 
of wave propagation. The classical waveform distortion process typical of intense acoustic waves in homogeneous 
media can be strongly altered when nonlinearly generated harmonics lie inside or close to band gaps. This allows the 
possibility of engineer a medium in order to get a particular waveform. Examples of this include the design of media 
with effective (e.g. cubic) nonlinearities, or extremely linear media (where distortion can be cancelled). The presented 
ideas open a way towards the control of acoustic wave propagation in nonlinear regime. 
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1. Introduction 

One of the most celebrated effects of wave propagation in periodic media is the appearance of forbidden propaga¬ 
tion regions in the energy spectrum of electrons, or band-gaps. Most of the physics of semiconductors, and therefore 
many electronic devices, are somehow based on this concept [j]]. In the late 80’s, these ideas where extended by 
Yablonovich and John 0 to light waves (electromagnetic waves in general) propagating in materials where the op¬ 
tical properties like the index of refraction were distributed periodically. These materials were named, by analogy 
with ordered atoms in crystalline matter, as photonic crystals. The typical scale of the periodicity is given by the 
wavelength. Actually, not only light but any wave propagating in a periodic medium may experience the same effects, 
and acoustic waves are not an exception. Sound wave propagation in periodic media has become very popular in the 
last 20 years in acoustics, after the introduction of the concept of sonic crystals 0 . Exploiting the analogies with 
other type of waves many interesting effects, as the mentioned forbidden propagation bands (band-gaps), but also 
focalization, self-collimation, negative refraction, and many others have been proposed. We consider in this paper 
the simplest case plane waves propagating in a ID structure, formed by a periodic alternation of layers with different 
properties. Depending on the context, such a structure has been named a multilayer, a superlattice (particularly in the 
context of semiconductors) or a ID phononic crystal (this include more exotic structures, as the granular crystal or 
lattice 0 ]). The huge majority of the studies considered so far have assumed a low-amplitude (linear) regime, neglect¬ 
ing the nonlinear response of the medium. Intense wave propagation in nonlinear periodic media, and in particular 
the case of sound waves, is almost unexplored. In this paper we present different examples of new phenomena related 
to sound wave propagation in ID periodic media, where each of the layer has a nonlinear quadratic elastic response. 
Nonlinear acoustical effects in such structure have been studied only in a few works. In 0 the harmonic genera¬ 
tion process is described in a fluid/fluid multilayered structure (water/glycerine), based in a nonlinear wave equation. 
Also, acoustic soli tons in solid layered nonlinear media have been presented in ]6[]. More recently, the complementary 
action of nonlinearity and periodicity has been considered in 0, where an asymmetric propagation device (acoustic 
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Figure 1: Layered acoustic system with two different layers and second harmonic generation scheme. Here the lattice period is a = a\ + « 2 - 


diode) was proposed. There, the nonlinearity and the periodicity act at different locations and its effect is considered 
separately. The effects discussed in this paper are the result of the interplay between nonlinearity and periodicity. 
Here we describe how the geometrical and acoustic parameters of the structure can be used to control the harmonic 
distortion processes in a multilayer. The conditions required to selectively act on the nonlinearly generated spectrum, 
and therefore manipulate the waveform in the desired way, are obtained and discussed. 

The theory presented here has been developed for fluid-fluid (scalar) structures, however the main conclusions are 
extendable to fluid-solid or to solid-solid multilayers, if particular conditions are given. Also, the main conclusions 
of this paper are independent on the regime of the waves (audible, ultrasound,...), and therefore on the size or scale of 
the structure. Specially interesting is the domain when ultrasound waves belong to the Terahertz regime, where these 
ideas may find a great potential. The progress in miniaturization and the technological development allows currently to 
create phononic multilayers at scales even in the nanometer range (each layer contains then a small number of atoms). 
This structures are usually made of semiconductors and are often used in particular applications as phononic mirrors 
to form phonon nanocavities @], or microcavities to obtain a strong optomechanical coupling 0 (for a revent survey, 
see (gtl). In a remarkable recent achievement, acoustic amplification was realized in doped GaAs/AlAs superlattices, 
where a SASER (Sound Amplification by the Stimulated Acoustic phonon Radiation) was demonstrated, in a device 
including a superlattice gain medium and GaAs/AlAs SLs acoustic mirrors 0 - 

The structure of the paper is as follows: In Sec. 2 we present the model for nonlinear propagation of acoustic 
waves in periodic media. The next Sec. 3 describes the process of harmonic generation in homogeneous media, 
and how it is modified by the presence of periodicity. In Sec. 4 the possibility of manipulating the spectrum of a 
propagating sound wave by tuning the parameters of the layered medium is discussed, showing examples of particular 
situation, as the case of a cubic-effective medium made out of quadratically nonlinear layers. Finally, Sec. 5 presents 
the conclusions. 


2. The model 

2.1. The medium and its dispersion relation 

We consider a periodic medium made of an arrangement of homogeneous fluid layers of thickness a\ and ao 
with different material properties. For the shake of simplicity only longitudinal waves under normal incidence are 
considered. A scheme of the medium is shown in Fig. [T] 

The propagation of small amplitude waves in an infinite periodic system is completely described by its dispersion 
relation, often known as band structure, that for ID systems as in Fig.|T|can be expressed analytically as 0 


cos (lea) = cos (k\a\) cos (£ 2 * 22 )-( —- H-I sin (k\a\) sin (^ai) 

2 \k 2 k\ 
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Figure 2: Dispersion relation of the two-layers system for layer proportion a = 0.5 and for different q/cj ratio. Left: real wavenumbers. Right: 
imaginary part of the complex wavenumber. 


also known as the Rytov formula, where k is the Bloch wave-number, a — a\ + ci 2 is the lattice period, and k, = w/c, 
is the local wavenumber, with c, the sound speed in the i layer. For a wave of frequency oj incident in a medium with 
known acoustical c, and geometrical a, parameters, the above equation results in a band structure of propagating and 
nonpropagating (bandgap) regions, as shown in Fig. [2] Thus, using Eq. (|7}, we can estimate the effect of periodicity 
on the different harmonics of the incident wave as they propagate through the multilayer, which is the main premise 
of this work. The ratio between layer thickness can be defined as a - a\/a, leading to aj = (1 — a)a. 

An example of dispersion relation plot is shown in Fig.[2]for normalized parameters a = 0.5 and for different sound 
speed ratios c\/c 2 - Increasing the impedance ratio between layers increases the reflected intensity in the trans-layer 
propagation, while the transmitted energy of the multiple internal reflections diminishes. As can be seen, due to these 
scattering processes band-gaps are progressively opened around the wavenumber k — nn/a with n = 1,2,.... Thus, 
the bandwidth of these band-gaps also increases when the impedance ratio grows. 

On the other hand, its imaginary part increases in amplitude with c\/c 2 , leading to shorter evanescent propagation 
inside the band-gap for high sound speed contrast layers, while remains zero (no attenuation) in the propagation 
band. We recall that the system is conservative: the physical interpretation of the complex wavenumber is not energy 
absorption, but back-reflection of the incident wave. Thus, at band-gap frequencies the waves penetrate only a short 
distance into the medium with a forward evanescent mode, and if the medium is perfectly periodic and lossless the 
energy is back-reflected (it behaves as a mirror). 

2.2. Nonlinear constitutive mode 

The nonlinear propagation of sound in the acoustic inhomogeneous media, and in particular in multi-layered 
media can be described by several models, with different levels of accuracy. Here, we use the equations of continuum 
mechanics for ideal fluids with space dependent parameters. These are the continuity equation for mass conservation 
0 


-f- + V • (pv) = 0. (2) 

at 

and the equation of motion that follows from conservation of momentum 

Dx 

p—+Vp = 0, (3) 

where p is the total density, v is the particle velocity vector over a Eulerian reference frame, p is the acoustic pressure, 
t is the time and D is the material derivative operator. 

For non homogeneous media, the ambient properties of the fluid in the absence of sound are space dependent, 
so the total density becomes p(t, x) = p'{t, x) + po(x), where po(x) is the spatially dependent ambient density and 
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p'(t, x) is the perturbation of the density or acoustic density, that is space and time dependent. Then, using the material 
derivative, Eq. (0 becomes 


P° ¥ + Vp = -p'—-(p'+ p 0 ) (v • V) v, (4) 

In this equation, the first two terms in the left-hand-side account for linear acoustic propagation, where the terms 
in the right-hand-side introduce nonlinearity in the Eulerian reference frame through momentum advection processes. 
On the other hand, we can expand Eq. © for nonhomogeneous media as 


— + poV • v + v ■ Vpo = -p'V ■ v - v ■ Vp'. (5) 

dt 

Here, the first two terms on the left-hand-side account for linear acoustic propagation, the third, also linear, ac¬ 
counts for the magnitude of the changes in the ambient layer properties. Note this term is space dependent but only 
changes at the interface between adjacent layers. For density matched layers, p,- = p, |, this terms vanishes. The terms 
on the right-hand-side are nonlinear and accounts for mass advection. 

Finally, a fluid thermodynamic state equation p = p(p, s ) is needed to close the system, with s the entropy. The 
local nonlinear medium response relating density and pressure variations, retaining up to second order terms, can be 
written as 


2 / B c o ,2 
P = C « P + 2A7 q P ’ 


( 6 ) 


where B/A(x) is the quadratic nonlinear parameter and co(x) is the sound speed, that can be also spatially dependent. 

In this system of equations, quadratic nonlinearity appears in the equation of motion © and in the continuity 
equation ©, in the momentum and mass advection terms respectively, and also in the equation of state, Eq. ©, re¬ 
lating pressure and density acoustic perturbations. We note that here we only take into account nonlinear processes 
through the layer’s bulk. The nonlinear effects at the boundary between adjacent sheets are neglected. These nonlinear 
boundary effects include cavitation processes, that in the case of fluids with very different compressibility can be very 
important. In the case of solid layers, other local nonlinear effects relative to boundaries, e.g. clapping phenomena be¬ 
tween surfaces, can lead to nonlinearities that are orders of magnitude in importance compared to the bulk cumulative 
nonlinearities. 


2.3. Second-order model 

For moderate amplitudes, the system of Eqs. ©0 can be simplified. For that aim, we use a perturbative method 
with same ordering scheme as in ill411 . where 0(e), 0(e 2 ) and 0(e 3 ) represents the terms of generic smallness param¬ 
eter e. The derivation of a second-order nonlinear wave equation requires the substitution of the linearized acoustic 
approximations (first order) into second order terms of Eq. 0©. This substitution procedure will give third order 
errors, so the final nonlinear wave equation will be a second order approximation of the full constitutive relations. 

These equations can be combined to form a single nonlinear wave equation valid for nonhomogeneous media up 
to second order approximation 


2 1 d 2 p 1 

^~P - — VpoVp = - 


dt 2 


Po 


P d 2 p 2 

Poc a q dt 2 


]_cP) 

c 2 0 dfi ) 


L + 0(s 3 ). 


(7) 


where we introduced the coefficient of nonlinearity /3 — 1 + ^ that accounts for material and mass advection quadratic 
nonlinearities. It is worth noting here that the second-order Lagrangian density vanish for plane progressive waves due 
to the first order relation p = ucopo that leads to X = 0. In this case, Eq. © simplifies to the well-known Westervelt 
equation for inhomogeneous media 


V 2 / 


1 d 2 p 1 


c 2 dt 2 


VpoVp 


P d 2 p 


, 2„2 


+ 0(E 3 ). 


( 8 ) 


Po PqCq dt 2 

In general, the Lagrangian density term can be discarded based on the distinction of cumulative and local non¬ 
linear effects. In this way, for progressive quasi-plane wave propagation in homogeneous media the nonlinear local 
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Figure 3: Harmonic generation in the layered medium at low frequencies (numerical results), and its comparison with analytical expressions 
(Fubini) for an homogeneous medium. 


effects become insignificant in comparison to the nonlinear cumulative effects, where in most practicals situations, 
beyond a distance of only few wavelengths away from the source local nonlinear effects can be neglected. However, 
local nonlinear effects can become significant in other complex situations including standing-wave fields and finite 
amplitude acoustic waveguides. Concerning the layered media, in this work we solve numerically the full constitutive 
relations, and the effect of the Lagrangian term is shown to be negligible under the conditions of our study. 

3. Harmonic generation in layered media 

We will study the response of the layered system for plane-harmonic wave excitation. Then, as sketched in Fig.Q] 
the source is placed at one boundary of the layered system, and the acoustic relevant magnitudes are calculated 
along space and time. As the wave propagates, cumulative nonlinear effects generate harmonics of the fundamental 
frequency, o>q, and due to the multiple scattering processes into the layers, local nonlinear effects also distorts the 
wave. However, the high dispersion of the layered system have a strong impact on the nonlinear harmonic generation. 
Dispersion modify the resonance conditions between fundamental and second harmonic wave, and also for other 
nonlinearly generated higher frequencies. In this way, nonlinear energy transfer efficiency from one component to 
another is modified in a wide variety of configurations, leading to the possibility of engineering and controlling the 
nonlinear wave processes by tuning the dispersion relation. 

Depending of the frequency of the input wave, different scenarios can be observed, as reported in the following 
subsections. 

3.1. Nondispersive (Fubini) regime 

We start studying the propagation in the layered system for harmonic excitation in the very low frequency regime, 
where we assume that ka <K 1 holds. As the Rytov’s Eq. 0 predicts, in the very low frequency regime the slope of 
the u>(k) curve is nearly constant. The dispersion of all the spectral components is negligible, and they all propagate at 
nearly the same velocity and are correspondingly phase-matched. Thus, in the absence of dispersion and attenuation 
process, the system of Eqs. ([2][3]) and <[6]> can be reduced for a harmonic-plane wave to a Burger’s evolution equation 
expressed in traveling coordinates with effective parameters, namely co, po and /?. An analytic solution of this equation 
in terms of the nth-harmonics of the fundamental wave of frequency a> and initial amplitude p o is known as the Fubini 
solution. 


p(cr, t) = po ) — J„ (ncr) sin (uwr) , (9) 

ncr 

n= 1 

where J„ is the Bessel function of order n, and cr = x/x s is the propagation coordinate, normalized to the shock 
formation distance, x s = 1 /[&. with the effective match number s = u/co and the effective wavenumber k = w/cb, 
that can be also found from Eq. ([]}. This celebrated solution is valid for cr < 1 (pre-shockregion). 
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Figure 4: Scheme of the phase miss-matching situation. The fundamental wave vector k\ at frequency oj generates a. forced wave 2k\ at frequency 
2to. The free wave that the system allows to propagate is k 2 , located in the dispersion relation curve. Due to dispersion, k 2 =£ 2k\, thus there exist a 
phase mismatch, A &2 between both waves and the generation is therefore asynchronous. 


Simulations were carried out using a full-wave constitutive relations solver. Thus, we shall define the normalized 
reference frequency as Qo — nco/a (located in the first band-gap). The source frequency was set to o) — 0.1 Q () . 

Figure [3] shows the analytical and numerical solutions for the low frequency limit of the layered system, where 
an excellent agreement is obtained between Fubini and numerical solutions in the pre-shock region, cr < 1 and for 
low excitation frequencies. As commented above, when the fundamental frequency is increased the higher harmonics 
fall in dispersive region of the frequency bands, and thus its wave speed is reduced. In this situation, phase matching 
conditions are no longer fulfilled and therefore the energy transfer from fundamental to higher harmonics is mod¬ 
ified. Thus, the Fubini solution can be only applied as an ideal solution for the low-frequency limit or as a good 
approximation for the first harmonics and for frequencies below o> < O.IOq. 

3.2. Dispersive regime 

For frequencies above the (idealized) homogeneous-Fubini regime, finite (weak and strong) dispersion effects are 
observed. The dispersive effects of the layered system deeply affects harmonic generation processes. 

As intense waves propagate through a quadratic nonlinear medium, their frequency components interact with each 
other and new frequencies arise at combination frequencies, including higher harmonics. The cumulative energy trans¬ 
fer from the interacting waves to the harmonics is dependent on the resonance conditions oii ±o >2 — k] ± k 2 = k 2 . 
Note these conditions express the laws of conservation of energy ( tuo ) and momentum (hk) in the quantum descrip¬ 
tion for the disintegration and merging of quanta Q. These conditions can be satisfied in a variety of situations. 
The most simple case is observed in nondispersive media and for collinear waves k, = ojfcQ. In this situation the 
resonance conditions are fulfilled all over the spectra and a large number of harmonics interacts synchronously: when 
there exist in the system a free wave with velocity /|k 3 1 that matches the excited (forced ) wave o>\ ± on/lki + k 2 |, 
th efree wave is excited in a resonant way. The resonant interaction leads therefore to synchronous (phase matched), 
cumulative energy transfer from the initial wave to the secondary wave fields. 

In the case of an initial monochromatic wave, the main wave generates its second harmonic. The resonant condi¬ 
tions in this situation read 2u>i - oji, 2k[ = k 2 , that holds true for nondispersive collinear waves, leading to the simple 
relation 2k(o>\) - k(2o>i). However, in the case of dispersive media this condition is, in general, not fulfilled and the 
forced and free waves interact asynchronously. Figure [4] shows such situation for a layered media with a fundamental 
wave in the first dispersion band. 

In order to study asynchronous second harmonic generation processes, we recall here for the lossless second- 
order wave equation Eq. ©, for one-dimensional propagation. This equation does not include dispersion by itself, 
dispersion arises from the solution of the linearized wave equation with the layered media boundary conditions, where 
the eigenvalue problem leads to the Rytov’s dispersion relation Eq. 0}. 
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In the following, we apply a perturbation method to obtain an approximate solution for the second harmonic field. 
We expand the pressure field as sum of contributions of different orders, i.e. p = p [l) + ep {2) + ■ ■ •, where s is the 
smallness perturbation parameter, which we identify with the acoustic Match number. Thus, p (l> is the first order 
(linear) solution of the problem and p t2) its the second order contribution. By substituting the expansion in the second 
order wave Eq. ©. assuming constant density!] and neglecting 0(s 3 ) terms we get a coupled set of equations that 
can be solved recursively. The solution of the first order equation corresponds to a monochromatic plane wave of 
frequency to 

p (l) - po sin(wf - k\x) (10) 

where k\ — k{to) is the wave vector associated with the primary frequency to, and po is the excitation pressure 
amplitude. Substitution of the first order solution into the equation obtained at the next order in the expansion, leads 
to an inhomogeneous equation for the second harmonic field: 


d 2 p (2) | d 2 p C) 

dx 2 Cg dt 2 


4/3u 2 P 2 0 

P0 C 0 


sin (2tot - 2k\x). 


( 11 ) 


The general solution of the this equation is the sum of the solution of the homogeneous equation (po = 0), and 
the particular solution of the inhomogeneous equation (p o ^ 0). Therefore the field for the second harmonic can be 
expressed as p (2) = pf 1 + pf\ where the corresponding waves for this two solutions are th e free, and forced waves 
respectively. Such homogeneous and particular solutions are: 


(2) 

Ph = 

pf ] \ x =o sin(2oV - k 2 x). 

(12) 

^3 

II 

n , o; v/ sin (2tO\t 2k lX ), 

(k 2 + 2k l )(k 2 - 2ki) 

03) 


where k 2 — k(2oo\) is the wavenumber of the free wave at second harmonic frequency, and the constant A = 
-A(5to 2 x p^l poCq. It is worth noting here that as long 2k\ + k 2 , the forced and free waves in dispersive media have 
different phase speed, i. e. the forced and free waves are phase mismatched as can be seen in the argument of the 
sin function in Eq. (1121113b . Imposing the boundary condition, that the second harmonic must be absent at x — 0, the 
second harmonic field can be expressed as 


P (2) = 


k 2 Ak 


A k \ 

—xl cos ( 2co\t - k' 2 x), 


(14) 


where the effective wave number is k' 2 - (k 2 + 2k\)/2 
second harmonic generation is defined as 


k 2 and the detuning parameter that describes the asynchronous 


Ak-ki- 2k\ = k(2to) - 2k(co). (15) 

Equation (IT-fl ) describes the well-known effect in second harmonic generation in dispersive media, that is the 
beatings in space of the second harmonic field when the resonant conditions are not fulfilled. Thus, as A k increases, 
the beating spatial period and also its maximum amplitude decreases. The position of the maximum of the beating, 
also called the coherence length, can be related to the second-harmonic phase-mismatching frequency as 


e |AA'| \k(2to) - 2k{to)\ 

This length corresponds to the half of the spatial period of the beating, where the maximum of the field is located. 
It can be expressed also for other higher harmonics simply as x c (n) = n/\Ak n \ = n/\k(nto) - nk(io )|. 

In the limiting case of A k —» 0, the second harmonic field is generated synchronously and accumulates with 
distance, so a linear growth is predicted. In this case, phase matching conditions are fulfilled and the free wave is 


1. We neglect the ambient density variations for the sake of simplicity. Dispersion arise also for sound speed variations, that are assumed to be 

implicit in the boundary conditions. 
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Figure 5: Second harmonic evolution for x c /x s = (1,1/4 and 1/8) obtained using Eq. < 1141 (continuous line), numerically (white circles), nondis- 
persive linear law of growth (dotted line) and Bessel-Fubini nondispersive solution (dashed). 



Figure 6: Evolution of the second harmonic field propagating in bang-gap for second harmonic frequencies (a) just above band-gap 2coq = 0.84£2o> 
(b) in the middle of the bandgap 2cjo = Ho. All results for a layered medium with a = 1/2 and c\/c 2 = 1/2. (c) Detuning of the second harmonic 
(continuous line) and imaginary part (dotted line) in function of the normalized frequency. 


excited synchronous to the forced wave. Note here that in the derivation of Eq. ( 1 1 4b only second order processes are 
taken into account and therefore, only second harmonic is predicted. This leads to overestimate the second harmonic 
field: as long no energy is transferred to third harmonic, second harmonic predicted by Eq. (IT4ll in the absence of 
dispersion grows indefinitely. The validity of this model can be explored expanding the Bessel functions of Fubini 
series near the source. A simple comparison between the full Fubini solution and linear second harmonic growth gives 
a reasonable approximation for distances cr < 0.5 or for second harmonic field values of p(2a>) < po/4-. 

Figure[5]shows three different simulations in the dispersive regime of the layered media where the wave amplitude 
and frequency has been selected to match x c /x s = 1,1/4 and 1/8. The higher beating spatial period waves corresponds 
to lower frequencies. The analytical solution for the second harmonic matches the full-wave numerical solution. 
However, differences can be observed in the second harmonic amplitude estimation for x c /x s = 1 (Fig. 0(a). This 
overestimation by the analytical solution can be related to the absence of energy transfer to higher harmonics, that is 
not considered by the perturbation solution but is included in the simulation and also in the Bessel-Fubini solution. 
Therefore, this model is specially suitable in situations where the third harmonic does not grow cumulative with 
distance. In the lossless layered media, this situations include frequencies that leads to very high-third harmonic 
detuning and also when the third harmonic falls in band gap. 

3.3. Second harmonic in band gap 

Waves with frequencies falling into the band-gap of the dispersion relation are evanescent due the non negligible 
imaginary part of its complex wave number. Thus, its amplitude decays exponentially with distance. If the nonlinearly 
generated second harmonic falls into a band-gap, its amplitude does not decay but reaches a constant value @]. Figure 
[6] shows this case for two different frequencies. The constant amplitude value of the second harmonic wave depends 
on the imaginary part of the wave vector. 

This effect can be understood in terms of the free and forced waves. If the second harmonic is evanescent (as 
follows from the dispersion relation), the wave will not accumulate with distance. The fundamental wave is “pumping” 
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energy to the second harmonic field at every point in space. Thus, the second harmonic field is generated locally and 
remains trapped inside the layered media. It reaches a constant level that depends on three main factors. In first place, 
the ’’pumping” rate, characterized by the fundamental wave amplitude and medium nonlinearity, or more strictly the 
ratio between the layer thickness and the shock distance a/x s . Secondly, it also strongly depends on the magnitude of 
the imaginary part of the complex wave number, i.e. the ratio between its characteristic exponential decay length and 
the shock distance in a layer. The characteristic decay length of the evanescent propagation is always shorter when the 
second harmonic is in the middle of the band-gap, leading to a weaker second harmonic field in this frequency region, 
as seen in Fig. [6] Finally, it depends also on the detuning of real part of the wave number, where for the first band-gap 
is minimum at the center. The first factor can be isolated and studied separately. However, the two last factors are 
linked through the specific dispersion relation of the medium. 

Figure [6jc) shows the detuning of the second harmonic and the imaginary part as a function of the frequency for 
a medium with a = 1/2 and c \ / 1:2 — 1/2, showing that at the middle of the band-gap these two factors have opposite 
effects: detuning is null (phase matching) when evanescent decay is nearly maximized, and viceversa. However, the 
magnitude of the effects can be very different. As the rate of the second harmonic generation (see the initial slope 
in Fig. 0 is independent on the detuning, and the evanescence implies that the wave decays after few layers, there 
not exist a practical compensation of the effects at the center of the band-gap. However, the situation is different for 
frequencies around the limits of the band-gap, where the coherence length is of the order of the exponential decay 
characteristic length. Thus, for frequencies just above bad-gap and for amplitudes with shock distance comparable 
to the evanescent characteristic decay length, the beatings can be also observed, as shown in Fig. [6] (a). Then, if 
frequency is increased the characteristic decay length becomes shorter than the shock wave distance and beatings 
cannot be observed, leading to to the characteristic constant second harmonic field shown in Fig. [6] 

3.4. Fundamental harmonic in band gap 

When the fundamental frequency of the wave lies within the band-gap, small amplitude waves propagate evanes- 
cently. Essentially, the same applies to finite amplitude harmonic waves. In general, if the shock distance is large 
compared to the characteristic decay length of the evanescent wave, the nonlinear effects have no chance to accumu¬ 
late and harmonic amplitude is negligible. Since the characteristic exponential decay is about few lattice sites, this 
means that the initial amplitude necessary to achieve nonlinear effects in this configuration is much higher than those 
in the preceding sections. Figure [7] shows the evolution of the first and second harmonic waves for a fundamental 
frequency at the Bragg frequency, 2«o = 1 Q lt . and with a frequency just above but into the band-gap, cjq = 0.87Oy 
for a layered media of a = 1/2 and c 1 / cp = 1/2. In the first case, the imaginary part of the wave vector is remarkable 
high and the waves decay fast after few lattice units. Due to this fast decay, the second harmonic interacts only over 
a short distance with the first, and its amplitude is very limited. After a few lattice units, the fundamental wave can 
be treated as a small-amplitude evanescent-wave. The second harmonic, that also falls in bandgap (but in the second 
band gap) also decays exponentially. 

On the other hand, if the fundamental frequency is set just above the band-gap, where the imaginary part of 
the wave-vector is smaller, the amplitude of the fundamental wave decays more slowly, penetrating deeper into the 
material. The interaction region with the second harmonic is larger, and nonlinear effects result in a more efficient 
generation of the second harmonic. Furthermore, as long the different (higher order) bandgaps in the layered media 
can have different bandwidth, in this configuration at loq = 0.87Oq second harmonic does not fall inside a bandgap. 
Therefore, the generated second harmonic wave at the beginning of the lattice propagates through the medium essen¬ 
tially without amplitude change. Due to the evanescence of the fundamental wave, there is only forced wave at the 
beginning of the medium. Therefore, although in this configuration waves are phase mismatched, beatings are not 
present: only the free wave propagates through the medium. 

4. Nonlinear acoustic field management 

4.1. Tuning nonlinearity with dispersion 

In the preceding sections we have explored the fundamental behavior of nonlinear waves generated inside the 
layered media. But also, medium parameters can be designed to provide specific conditions. The material parameters 
can be tuned to get coherence at one frequency of interest, e.g. at one of the harmonics of the fundamental wave, or 
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Figure 7: (a) Evolution of the fundamental harmonic wave field with its fundamental frequency falling just above into band-gap, cjq = 0.87Qo> 
(continuous line), and in the middle of the bad-gap, 2cjq = lf^o (dotted line), (b) Corresponding second harmonic field, where for ojq = lQo 
(dotted line) second harmonic frequency falls in the 2nd band band-gap while for ojq = 0.87Qo> (continuous line) lies into a propagating band. 


to get detuning or evanescent propagation at other specific harmonics. Using these mechanisms the layered medium 
can be used to provide a balance of the harmonic amplitudes, or to obtain specific nonlinear waveforms, providing a 
control of the nonlinear process inside the medium. 

In the design of a system for this purpose, the coherence length is a useful control parameter. For this aim, the 
analytic Eq. (QJ is used, which is shown to provide an excellent framework to tune the layered parameters to obtain the 
desired balance between detuning, evanescent propagation, synchronous generation and, at the same time, it allows 
to find those conditions for a specific phase/group speed. Figure |9ja,b) shows an example of a dispersion relation, the 
coherence length for the second and third harmonic. The resulting harmonic amplitudes when phase matching of all 
harmonics is achieved is shown in Figure[8] This happens for a set of frequencies ojq = (0,1.75,2.333, ...)/Q. On the 
other hand, there also exist frequencies at which there exist coherence for the second but a non-negligible detuning is 
observed for the third. The opposite effect can be also obtained, where coherence is achieved for the third harmonic 
but second harmonic presents strong dispersion. Finally, other interesting regions are those where second harmonic 
component is almost phase matched and for the same frequency third harmonic falls into a band-gap. 

In the following subsections, we propose and analyze different configurations of the layered medium with specific 
balance between detuning, evanescent propagation and synchronous generation. 

4.2. Enhanced second harmonic generation 

One can expect that second harmonic generation is maximized in homogeneous nondispersive media. However, 
in nondispersive media coherence is achieved not only at second harmonic frequency, but also in the higher spectral 
components. As a result, energy is transferred from second harmonic field to higher spectral components and therefore 
second harmonic field does not grow indefinitely. Moreover, shock waves are formed and nonlinear absorption reduces 
wave intensity for cr > n/2 even in lossless media 0. 



Figure 8: Harmonic distribution for the frequency cj o = 1.75/H. Coherence is recovered for at least the lowest spectral components. Blackstock 
solution (dotted lines). 
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Figure 9: (a) Dispersion relation for a layered media of c\/c 2 = 1.33 and a = 1/2. (b) Coherence length for second (red) and third (blue) harmonics 
as a function of the fundamental frequency. Phase matched frequencies are those with x c —> oo, while asynchronous generation is predicted for 
x c —> 0. Frequencies at which the fundamental frequency is in band-gap are are marked in gray regions, while band-gap regions for second and 
third harmonic are marked in dashed lines.(c) Harmonic distribution for cj o = 1.668£2o where a coherence is achieved for second harmonic while 
the frequency of the third harmonic falls into the bad-gap. (d-g) (continuous lines) Waveforms at different distances for cj o = 1.668Qo. At cr = 3.3 
second harmonic generation field is maximize and can be seen the period doubling in the waveform. Then, at cr = 7.8 due to second harmonic 
detuning nearly sinusoidal wave is recovered. Analytic Fubini-Blackstock solution for the harmonics (red dotted lines) are plotted for comparison. 


The dispersion of the layered system can be used to modify this situation by including phase mismatches that alter 
the higher harmonic cascade processes, while maintaining coherence for the second harmonic. Figure [ 8 ] shows an 
example of a dispersion relation where for om> = I. 668 Q 0 it can be observed that there exist a reasonable coherence 
for the second harmonic ( x c /a = 1000), while the third harmonic falls in a band-gap. Figure [9] shows the harmonic 
distribution in this situation. Here, energy is transferred to second harmonic field that grows almost linearly for cr <2. 
On the other hand, the energy transferred from second to third harmonic is not cumulative and its amplitude does 
not grow with distance. Third harmonic experiment evanescent propagation due to the imaginary part of the complex 
wave-vector at this frequency. A constant field, as studied in Sec. 13.31 is obtained for the third harmonic. 

The total amount of the second harmonic amplitude in nondispersive media is /? 2 lmax ~ 0.36po, while in the 
example of Fig. [9] a maximum second harmonic amplitude of piLax ~ 0.75/?o is predicted. As can be shown the 
decreasing of the first harmonic follows the analytic nondispersive Blackstock solution for cr ;$ 3. Thus, in this 
regime all the energy of the first harmonic is being transferred to the second harmonic field. However, due to finite 
detuning of the second harmonic a long spatial beating is produced, with normalized period 8 cr, and energy is returned 
back to the first harmonic component. 

It is worth noting here that at distance cr « 3 sawtooth profile is observed in the nondispersive media. In contrast, 
only second and first harmonic have remarkable amplitude into the layered media. Waveforms are shown in Fig. [9j d- 
g). Near the source, where the amplitude of higher harmonics in not relevant the nondispersive waveform (in red 
dotted) is well approximated by the fundamental and its second harmonic of the layered medium. However, due to 
the evanescent propagation of the third harmonic for longer distances the nonlinear solution of the layered medium is 
mainly composed by the fundamental and its second harmonic. The maximum second harmonic in this configuration 
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Figure 10: (a) Dispersion relation for a layered media of c\/c 2 = 1/3 and a = 0.3. (b) Coherence lengths for second (red) and third (blue) 
harmonics, (c) Harmonic distribution for ojq = 0.529 /Q. where a coherence is achieved for second harmonic while the frequency of the third 
harmonic falls into the bad-gap. Bottom: (continuous lines) Waveforms at different distances for co o = 0.529 /Q.. At cr = 3.3 second harmonic 
generation field is maximize and can be seen the period doubling in the waveform. Then, at cr = 7.8 due to second harmonic detuning nearly 
sinusoidal wave is recovered. Analytic Fubini-Blackstock solution for the harmonics (red dotted lines) are plotted for comparison. 


is observed at cr = 3.3, as it can be appreciated in the waveforms of Fig. [9] the period doubling. Moreover, due to 
finite detuning of the second harmonic the process is not cumulative for all distances and at cr = 7.8 the energy is 
restored in the first harmonic again and a sinusoidal wave is obtained. Note that not all the energy is restored to the 
first harmonic in Fig.[9]at cr = 7.8, leading to a sinusoidal wave of different amplitude as can be observed in Fig. [9] 
The energy loss is mainly due to the artificial (numerical) viscosity necessary to nonlinear convergence |3- For these 
simulations the total distance is 1200 lattice sites and therefore the effects of attenuation are not negligible. However, 
the main nonlinear effects related to strong lattice dispersion still appreciated. An analogous effect has been also 
studied 01311 where instead of dispersion, selective absorption at specific frequencies is used to modify and enhance 
harmonic generation. 


4.3. Enhanced third harmonic generation 

In the first band (co < Qo), coherence is always lower for the third harmonic than for the second. However, in 
the superior bands the layered medium parameters can be tuned to obtain higher coherence for the third than for the 
second harmonic. Essentially we follow same ideas on the preceding section but for the third harmonic. In this case, 
the lattice is designed forcing the second harmonic to fall in bandgap. A the same time, perfect coherence can be 
found for the third harmonic at u> — 1 .4Qq. This situation is illustrated on Fig. [TO] around oj = 1.4Q () . In this case, the 
dispersion relation was obtained for a layered medium with parameters a = 0.3 and cijc\ = 1 /3. 

In this situation, as FigITOlshows. the second harmonic wave attains a constant value of about 0.04/?o- As discussed 
in Sec 13.31 this constant field does not grow with distance and is related to the evanescent solution of the free wave 
and the local nonlinear “pumping”. On the other hand, due to the coherence of the third harmonic, all the energy 
transferred form second to third is accumulated with distance. Therefore, near the source the rate of energy transfer 
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from second to third harmonic is constant. Thus, third harmonic start to grow almost linearly with distance, opposite 
to quadratically in homogeneous media. 

Numerical simulations also show fourth and fifth harmonics grow (not shown in Fig.lTOl). but only fifth harmonic 
harmonic reach a remarkable amplitude, growing near the source almost quadratically with distance. Therefore, the 
entire system behaves as an artificially cubic-like nonlinear medium formed by quadratic nonlinear layers. 

The corresponding waveforms measured at cr = (0.5,2,8,12) are shown in Fig. 1 1 Ot d-g). For <j = 0.5 and 2, it can 
be observed how the wave steepens with the characteristic shape of cubic nonlinearity. No shock waves are formed as 
long as strong dispersion is present for high frequency harmonics. It is worth noting here a remarkable fact: it steepens 
in the positive time axis direction (to the right in the figure), opposite than the quadratic nonlinearity plotted in red 
dotted as a reference. This effect, i.e. the steepening on the opposite side of the propagation direction, is characteristic 
of materials with negative parameter of nonlinearity. Therefore, the effective nonlinear behavior observed by the 
simulations in this conditions can be described as negative-cubic-like nonlinearity. 


5. Conclusions 

The interplay of dispersion and nonlinearity in multilayered periodic media, such as one-dimensional phononic 
crystals or superlattices is shown to have a strong impact on the acoustic waves propagating through the structure. 
Nonlinearly generated harmonics propagating at different velocities are phase-mismatched, modifying the transfer of 
energy between the different harmonics, and therefore the waveform itself. Shock formation, typical of nonlinear 
homogeneous media, is in this way avoided. We propose a model and some particular solutions to study this problem, 
and report examples of configurations that result in an effective control of the spectrum of nonlinear acoustic waves 
by tuning the dispersion relation of the medium. Selective enhancement of second or third harmonic is demonstrated, 
leading in some cases to situations where the structure behaves with an effective nonlinearity different from that of its 
constitutive elements. 

The work was supported by Spanish Ministry of Economy and Innovation and European Union FEDER through 
project FIS2011-29731-C02-02. A. Mehrem acknowledges Generali tat Valenciana the support from Santiago Grisolia 
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